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Abstract 

—  For  systems  of  reaction-diffusion  equations  with 
Neumann  boundary  conditions,  it  is  shown  that  the  solutions 
are  asymptotic  to  the  solutions  of  an  ordinary  differential 
equation  if  the  diffusivity  is  large.  The  methods  apply 
also  to  reaction-diffusion  systems  with  time  delays. 
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Introduction  and  Statement  of  Results 


Many  models  of  chemical,  biological  and  ecological 
problems  involve  systems  of  reaction-diffusion  equations 
of  the  form 


(1.1) 


3u/3t  =  DAu  +  f (u)  in  (2 


3u/3n  =0  in  3(2 


where  (2  is  a  bounded  open  set  in  ]R  with  3(2  smooth, 
N 

u  £  ]R  ,  D  =  diag(dlf . . . ,dN) ,  where  each  dj  >  0  is  a  con 
N  N  2 

stant  and  f :  ]R  -*■  3R  is  a  C  -function.  Other  types 
of  boundary  conditions  may  also  occur. 

In  recent  years,  there  have  been  many  investigations 
devoted  to  the  study  of  stable  patterns  for  (1.1) y  that 
is,  stable  solutions  which  are  spatially  dependent  (see, 
for  example,  [1-8]).  For  the  understanding  of  how  stable 
patterns  are  created,  it  is  obviously  of  interest  to  char¬ 
acterize  those  situations  for  which  stable  patterns  do 
not  exist  and,  even  more  particularly,  those  systems  for 
which  the  flow  is  essentially  determined  by  the  ordinary 
differential  equation  (ODE) 


(1.2) 


du/dt  =  f (u) . 


This  latter  problem  has  been  investigated  by  Conway, 

Hoff  and  Smoller  in  [9]  under  the  assumption  that  there  is 

N 

an  invariant  region  £  for  (1.1)  .  A  set  £  c  ]R  is  an 


invariant  region  for  (1.1)  if,  for  any  initial  data  Uq 
with  uQ(x)  e  £  fo r  all  x  £  ft,  one  has  the  solution  u(t,x) 
through  uQ  at  t  =  0  also  in  E  for  all  x  £  ft.  In 
[10],  [2],  it  is  shown  that  an  invariant  region  E  for 
(1.1)  must  be  a  rectangle  in  3R  if  all  the  diffusion 
coefficients  dj  are  distinct.  If  the  diffusion  coefficients 
are  equal ,  then  E  can  be  any  convex  set  which  is  positively 
invariant  for  (1.2). 

Suppose  E  is  an  invariant  region  for  (1.1),  M  = 
sup{ | 3f (u) /9u j  :  u  €  E},  -  A  is  the  first  nonzero  eigen¬ 
value  of  the  Laplacian  with  homogeneous  Neumann  boundary 
conditions  on  ft,  d  =  min (d^, . . . , dN) ,  and  o  =  dA  -  m.  In 
[9] ,  it  was  shown  that  o  >  0  implies  the  solutions  of 
(1.1)  with  initial  data  in  E  approach  a  solution  of  the 
ODE  (1.2)  as  t  -*■  00 .  Since  A  is  inversely  proportional 
to  the  squared  diameter  of  ft,  the  hypothesis  a  >  0  says 
that  diffusion  on  the  domain  is  fast  relative  to  the  reaction 
term  f.  In  fact,  the  estimates  in  [9]  show  that  the  spatial 
inhomogeneities  are  quickly  damped  out  if  o  is  very  large. 

As  remarked  earlier,  the  hypothesis  that  E  be  an  in¬ 
variant  region  for  (1.1)  severely  limits  the  types  of  equa¬ 
tions  that  can  be  considered.  In  fact,  as  pointed  out  by 
Smoller  [2;p.212],  the  property  that  E  is  an  invariant 
region  for  (1.1)  is  not  continuous  with  respect  to  the  d^ . 

In  fact,  if  E  is  a  convex  positively  invariant  set  for 


Jt  Ijt  V*|, 


-  3  - 

the  ODE  (1.2)  and  each  d^  =  1,  j  =  1,2,...,N,  then  E  is 
invariant  for  (1.1).  However,  if  E  is  not  a  rectangle, 
then  it  is  not  invariant  for  (1.1)  unless  the  d..  remain 
equal . 

It  is  the  purpose  of  this  paper  to  begin  an  investiga¬ 
tion  of  the  behavior  of  the  solutions  of  (1.1)  when  the 

constant  dX  is  large  and  the  equation  (1.1)  may  not  have 

N 

an  invariant  region  in  H  .  The  methods  will  use  proper¬ 
ties  of  the  flow  defined  by  (1.1)  in  function  space  and 
will  be  applicable  to  situations  where  the  equation  with 
no  diffusion  is  a  retarded  functional  differential  equation 
or  a  differential  equation  with  delays.  In  this  latter 
case,  it  is  almost  impossible  to  have  an  equation  with  an 
invariant  region. 

We  now  describe  the  results  in  some  detail.  Let  X  = 

L2(^,1RN  )  ,  D(A)  =  {(j>€  W2,2(fl,3RN  )  :3<J>/3n  =  0  on  3^}, 

A  =  -A  :  D (A)  -*■  X.  In  the  usual  way,  one  defines  the 

fractional  power  spaces  xa  using  the  operator  A.  If 

3/4  <  a  <  1,  then  it  is  known  that  XQ  c:  W^'^(n,]RN  )  fl 
oo  n 

L  (ft,nR  )  with  continuous  inclusion.  One  can  then  show 
(see,  for  example,  Henry  [10,p.75])  that,  for  any  uQ  €  Xa, 
a  >  3/4,  there  is  a  unique  solution  u(t,*,uQ)  €  Xa  of 
(1.1)  through  u^  at  t  =  0  which  is  continuous  in  t,Ug. 
A  set  ja/  <=  Xa  is  a  compact  attractor  for  (1.1)  if 


ja/  is  compact,  invariant  and  there  is  a  neighborhood  U  of 
ja/  such  that  the  w-limit  set  of  U  is  _sa/ .  By  the  co-limit 


'-'s777Jr777?.''?rr 


■T  r_ 


-  “  -  -  ,-v  ii  Vi  !TV  T-  '-■»■ 
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set  w(U)  of  U,  we  mean 

ui(U)  =  H  Cl  U  u(t,*,U)  . 

T>0  t>T 

N 

In  a  similar  way  one  can  define  a  compact  attractor  in  3R 
for  the  ODE  (1.2) . 

If  is  a  compact  attractor  for  the  ODE  (1.2),  then 

it  can  be  considered  as  a  subset  of  the  constant  functions 

.  CL 

m  X  and  it  will  be  a  compact  invariant  set  for  the  PDE 
(1.1).  However,  it  need  not  be  an  attractor  without  some 
conditions  on  D  and  ft.  If  d  =  mintd^ ,  j=l,2,...,N)  and 
-A  is  the  first  nonzero  eigenvalue  of  A  with  Neumann 
conditions,  then  the  main  result  of  the  paper  is  the  follow¬ 
ing 


Theorem  1.1.  Suppose  is  a  compact  attractor  for  the 

ODE  (1.2).  Then  there  exists  a  6  >  0  such  that  con- 

.  ot 

Sidered  as  a  subset  of  the  constant  functions  in  X  , 

a  >  3/4,  is  a  compact  attractor  for  the  PDE  (1.1)  if 

dA  >  6 .  More  precisely,  there  is  a  neighborhood  V  of 

.  ct 

jq/  in  X  and  constants  K  >  0,  c  >  0  such  that,  for  any 
Uq  €  V,  the  solution  u(t,*,Ug)  of  (1.1)  through  u^  at 
t  =  0  satisfies 


| u ( t ,  •  , Ug )  -  u(t)  |  a  <  Ke  Ct,  t  >_  0 

X 


where  u(t)  =  I  ft 


,-l 


u(t,x,Ug)dx  and  u(t)  satisfies 


the  equation 


*£*  «;  ^ |k  r  , 
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du(t)/dt  =  f  (u ( t) )  +  g(t,uQ) 


where  |g(t,uQ)  |  _<  Ke  ct,  t  ^  0. 

The  technical  part  of  the  theorem  states  that  the  solu¬ 
tion  u(t,*,uQ)  approaches  its  average  value  u(t)  ex¬ 
ponentially  as  t  -*  00  in  the  space  Xa.  Since  a  >  3/4, 
this  implies,  in  particular,  that  the  solution  approaches 


—  00  N 

u(t)  exponentially  as  t  ■*  °°  in  L  (fi,]R  )  .  The  conclu¬ 
sions  in  Theorem  1.1  are  the  same  as  the  ones  in  Conway, 

Hoff  and  Smoller  [9]  mentioned  above. 

The  proof  of  the  theorem  is  given  in  Section  2  and 
uses  elementary  properties  of  Liapunov  functions  for  ODE's 
and  a  special  decomposition  of  (1.1).  In  Section  3,  we 
consider  invariant  regions  and  show  how  the  method  presented 
here  gives  the  same  qualitative  results  as  in  [9] ,  but 


the  rates  of  decay  are  not  as  sharp.  Generalizations  to 
functional  differential  equations  are  given  in  Section  4. 


2 .  Proof  of  Theorem  1.1. 

N  N 

If  M  c  ]R  is  a  given  set  and  x  £  I  is  given, 

we  let  d(x,M)  denote  the  distance  from  x  to  M.  Suppose 

^  is  a  compact  attractor  for  the  ODE  (1.2).  From  the  book 

of  Yoshizawa  [11, p. Ill],  there  is  a  neighborhood  U  of 

.0/  and  a  Lipschitz  continuous  function  V:  U  ■*  ]R  such 

that,  for  any  x  €  U, 


v  ■i-.*  ? - •*  A 


-  6  - 


(i)  V  (x)  =0  if  x£i 


(ii)  a(d(x,J^))  £  V(x)  £  b(d(x,_Q^))  where  a(r)  is 
continuous,  nondecreasing,  a(r)  >  0  if  r  >  0,  and  b(r) 
is  continuous,  b(0)  =  0. 

(iii)  V(1  (x)  £  -  V(x)  where 

2)(x)  =  lim^o  h  1  [V  (u  (h,x) )  -V(x)  ]  with  u(t,x) 

being  the  solution  of  (1.2)  through  x  at  t  =  0. 

In  the  following,  for  any  c  >  0,  we  let  V  =  {x  €  U: 

c 

V(x)  <  c>,  vc  ~  CAVc‘  From  property  (ii)  above,  Vc  is 

compact  for  any  c  >  0. 

a 

Let  W  c  x  be  the  linear  subspace  consisting  of  the 

constant  functions,  Xa  =  W  ©  W^,  u  =  v  +  w  where  v  6  W, 

w  €  W1, 
a 

(2.1)  v  =  |fi|  1  u(x)dx,  f  w(x)dx  =  0 

Jo  Jo 


We  can  identify  W  with  IR  and  therefore  will  consider 

N 

v  as  an  element  of  W  as  well  as  a  vector  in  2R  . 

Suppose  u(t,*)  is  a  solution  of  (1.1)  and  let 

u(t,*)  =  v ( t )  +  w ( t , • ) ,  v ( t )  €  W,  w(t, • )  €  W1.  Then 

a 


(2.2) 


dv/dt  =  P(v,w) 


3w/3t  =  DAw  +  Q(v,w) 


where 


(2.3) 


P  (v,  <j>)  =  Ini'1!  f  (v  +  4>(x))dx 

n 

Q(v,4»)  (x)  =  f(v  +  <J>(x))  -  Inf'1  f  ( v  +  4>{y))dy. 

o 


Observe  that 


(2.4)  P(v,0)  =  f (v) ,  Q(v,0)  =  0  . 


Also,  for  €  W^,  we  have 

[Q^  (v,4>)<i>]  (x)  =  f '  (v+4>(x)  )^(x)-|Q|-1J  f  ’  (v+<My) )  My)dy 
=  [f*  (v+*(x))-f'  (v)  ]  ^  (x)  +  f'(v)iMx) 


-  \n\~1\  [f '  (v+*(y))  -  f'(v)]iMy)dy 

since  ^  =  0.  Since  a  >  3/4,  <J>  €  Xa  implies  <j»  €  L  (fl, 

and  there  is  a  constant  k  such  that  |  <p  |  ^  <  k  1 4>  I  n-  Let 

_  LX 

M  =  sup{|f'(v) | :  v  €  V  }.  For  any  c.  <  c,  there  is  a 

C  C  X 

6  >  0  such  that  |  <J>  |  <  <5,  v  €  V  implies  v  +  4>(x)  e  V 

Xa  2  Cj- 

for  x  €  SI.  Since  f  is  a  C  -f unction,  it  follows  that 

there  is  a  constant  N  such  that 

c 


|f  •  ( v+  4*  ( x )  )  -  f  •  (V)  I  <  N  U|  <  N  k|  <*.  I 

C  L  c  x“ 

v  €  V  ,  |<j)|  <6.  Therefore, 

C1  x“ 

|  IQ^(V,*M>]  (x)  -  f '  (v)  (x)  |  <  2k2Nc  |  <j)  |  aU|  a 

X  X 


which  implies 


(2.5)  |Q^(v,*)<»  -  f '  (v)^|x  <  2k2Nc  |  «  | 15 1  <f>  |  Ji|>|  a 

X  X 

Relation  (2.51  and  the  triangle  inequality  imply  that 

I Q ( v , <j> )  |  <  e|<t>| 

x“ 

e  =  m  k  +  2k2N  |n|1//2s 

P  P  1  ' 


(2.6) 


9 


in  the  set  v  €  V 


d>  <  6 . 
xa 


■  vh 


Let  T  (t)  be  the  semigroup  on  W^,  0  <_  a  <  1,  generated 
by  the  equation 


3w/3t  =  DAw 
3w/3n  =  0 


in  ft 

on  3ft , 


Now,  fix  a  >  3/4.  There  is  a  constant  k1  >  0  such  that 


(2.7) 


| T ( t) w |  a  <  kjeT^^w)  a,  t  >  0,  w  6  w£ 
X  X 

|T (t) w|  a  <  k^~^Xtt~a |w|x»  t  >  0,  w  €  W1 
X 


We  can  now  use  the  variation  of  constants  formula  to  rewrite 
the  initial  value  problem  for  (2.2)  as 


dv/dt  =  f(v)  +  [P(v,w)  -  f(v)),  v(0)  «  vQ, 

(2.8)  t 

w ( t )  =  T(t)wn  +  |  T  (t-s)  [Q(v(s)  ,w(s) )]  ds. 

J0 

For  any  c.  <  c,  choose  A  as  before  so  that  |$|  <  S, 

_  X 

v  €  V  implies  v  +  4>(x)  €  V  .  Choose  a  constant  o  , 
cl  c 

0  <  a  <1,  let  k2  be  the  Lipschitz  constant  for  V  on 
V  ,  L  =  [  s"ae"U  0y/,X^sds  and  let  V(v)  >.  p  >  0  for 

c  _  Jo 

v  €  V  V  . 


•'.‘A 

.'.V 


With  the  constants  as  in  (2.6),  choose  dX  >  0  and 
so  that 


n  <  6 


(2.9) 


? d§fk1L6 (dX) 0-1  <  1 
p  -  kjM^n  >  0 


Let  v(t),w(t)  satisfy  (2.8).  If  v(s) 

|w(s)  |  <  n  for  0  £  s  <_  t,  then  relation  (2 

X 

V(v(t))  <  -  v(v(t))  +  k2|P(v(t)  ,w(t))  - 

<  -  V(v(t) )  +  k2Mck|w(t) |  o 

X 

<  -  v (v (t) )  +  k,M  kn 

—  2C 

z(t)  <  k  e"d(X-&)  tz  (0)  +  k.e  ft(t-s)"“ 

J0 

where  z(t)  =  |w(t) |  expdat 

Xa 

If  y (t)  =  sup{ z ( s ) ,  0  <  s  <  t},  then 


z(t)  <  k1e“d^X  °^tz(0)  +  Cy(t). 

This  implies  that  y(t)  <  (1-?)  -1z  (0)  . 

that 


|w(t)  | 


X 


a 


l-S 


“dor  t 
e  -  |w 


0  Xa 


.6)  implies  that 
f (v(t)) | 

e-d(X-o)(t-s)z(s)ds 

This  implies 


Thus,  we  see  that 


10 


(2.10) 


V (v  (t) )  <  -  V (v  (t) )  +  k,M  kn 

~  2  C 


|w(t)  |  <  — ^ —  e 

Y«  1-C 


-dot  i 


w, 


0  xa 


if  it  is  assumed  that  v(s)  €  V  ,  |w(s) |  <  n. 

C1  Xa 

Now  choose  6^  so  that  k^d-C)  <  n.  If 

vQ  €  Vc  ,  | wQ J  Q  ^  then  relations  (2.9),  (2.10)  imply 

that  v ( t )  €  V  ,  | w ( t) |  <  n  for  all  t  >  0. 

C1  X 

Relation  (2.10)  also  implies  that  |w(t) |  approaches 

Xa 

zero  exponentially  as  t  -*■  Thus,  the  w-limit  set  of 

every  solution  of  (2.2)  with  initial  value  v_  €  v  and 

u 

|wQ|  a  <  6^  must  have  w  *  0?  that  is,  the  w-limit  set 
X 

of  any  solution  of  (1.1)  with  uQ  *  vQ  +  wQ  satisfying 

the  above  condition  must  lie  in  the  set  W  n  v  .  Further- 

C1 

more 


dv ( t ) / dt  =  f (v(t) )  +  [P(v(t) ,w(t))  -  f (v(t) ) ] 


and  the  second  term  approaches  zero  exponentially  as  t  +  00 . 
Therefore,  the  limit  set  of  v  must  be  a  union  of  invar¬ 
iant  sets  of  dv/dt  «  f(v)  which  belong  to  V  (see 

C1 

Yoshizawa  [11]).  However,  all  such  invariant  sets  must 
belong  to  jaf.  This  completes  the  proof  of  the  theorem. 


3.  The  Case  of  an  Invariant  Region. 

If  equation  (1.1)  has  an  invariant  region  E  in 
®N  ,  and  v  €  E,  4>  €  w£,  v  +  4>(x)  €  E  for  x  €  fl,  then 


Conway,  Hoff  and  Smoller  [9]  prove  a  more  global  version  of 


Theorem  1.1.  In  fact,  assuming  the  initial  data  uQ  € 

ck'^(fi,IRN),  Ug(x)  e  £  for  x  €  fi,  they  prove  that  the 

solution  u(t,x)  through  Ug  satisfies  |Vxu(t,*)|  2' 

L 

(u(t,*)  -  u(t)  1  ^  0  as  t  ■*  “  exponentially  with 

L  2 

exponent  a.  To  obtain  the  L  estimate  on  V^u  and  an 
2  .  — 

L  estimate  on  u(t,*)  -  u{t)  is  not  difficult  using  simple 

integration  by  parts  and  the  variational  characterization 

00 

of  the  first  eigenvalue  of  A.  The  L  estimate  uses  more 
sophisticated  properties  of  parabolic  equations  and  the 
readfer  is  referred  to  [9]  for  details. 

The  method  used  in  the  previous  section  can  be 
applied  to  this  case  but  the  exponential  estimate  is  not 
as  precise.  We  briefly  indicate  how  this  can  be  accom¬ 
plished.  If  uQ  €  Xa,  a  >  3/4  and  uQ(x)  e  £  for  x  €  ft, 
then  the  solution  u(t,x)  through  Ug  remains  in  I  for 
all  t  >  0.  If  4>€W^,v  +  <j>(x)  G  E  for  x  6  ft,  then 
one  easily  shows  that  iQ^v,^)  |  <  2|fi|^kM  where  k  is  such 

that  |  4>  |  <  k  |  4>  |  and  M  =  sup{  |  f  '  ( v)  |  ,v  e  E  } .  Using 

L  “  Xa 

the  variation  of  constants  formula 


w(t)  =  e 


bwn  +  [  eDA(t~S^Q(v(s) ,w(s) )ds 
u  i  0 


and  letting  k^  =  2k^  |  fi |  ,  one  obtains 


l«(t)  Ix0  <  3Cie_<^Xt  lw0  l^ct 


+  2k.kMf  e"dX(t's) (t-s) ~a Iw (s) |  ds 


Let  z (t) 
Then 


z(t)  <  k^CO)  +  2k3kM|  e"(dX"B)  (t_s)  (t-s)“az(s)ds 

f0D 

<  k1z(0)  +  2k3kMj  e“(dX_e)ss"ads  y(t) 

It  is  clear  that  one  can  choose  g,d  so  that  the  coeffic 
ient  0  of  y(t)  is  <1.  Then, 

|w(t)  |  a  <  k^i-erV^lwQ^,  t  >  0  . 

X 

This  estimate  implies  that  |w(t)  |  and  so  |w(t)  | 

X  L 

exponentially  as  t  00 .  This  also  implies  that 
P  (v(t)  ,w(t) )  -  f(v(t))  -*•  0  exponentially  and  we  obtain 
all  of  the  conclusions  of  Theorem  1.1.  Of  course,  the 
exponential  estimate  is  not  a  good  one. 


4.  Functional  Differential  Equations. 

N 

Suppose  r  >  0  is  a  given  constant,  C([-r,0],IR  ), 
f  :C(  (-r,0]  ,IRN  )  -*■  ]RN  is  a  given  C2-f unction  and  consider 
the  equation 


(4.1) 


3u(t,x)/3t  =  DAu(t,x)  +  f(ut(*,x))  in 

9u(t,x)/3n  =  0  on  3fi 


where  u.  (0,x)  =  u{t+0,x),  e  e  [-r,0],  x  €  fi. 


As  before,  let  X  =  L  (Q,]R  ),  define  the  spaces  X 
and  choose  3/4  <  a  <  1.  For  any  4>  €  C(  [-r,0]  ,Xa)  with 
3<j>  (0  ,x)/3n  =0,  -r  <_  0  <_  0,  one  can  define  a  solution 
u(t,x,<j>)  of  (4.1)  on  an  interval  -r  £  t  <  a,  ot  >  0, 
with  u(0,x,<J>)  =  <j>{0,x),  "r  £  0  <_  0.  This  solution 
u  (t,  • ,  <J>)  will  be  in  xa,  it  will  be  continuous  in  t,4>  and 
continuously  differentiable  in  $ .  The  proof  of  this  fact 
follows  along  the  lines  of  Henry  [10]  or  Travis  and  Webb 
[11]  . 

If  we  let  T(t,<j>)  €  C([-r,0],Xa)  be  defined  by 
T (t,  <*>)  (0)  =  u  (t+0  ,*,<)>)  ,  -r  £  0  <_  0,  then  T  (t ,  • )  defines 
a  local  semiflow  on  C ( [-r, 0] ,xa) .  Positive,  negative  and 
complete  orbits  are  defined  in  the  usual  way  as  well  as 
a-limit  sets,  w-limit  sets,  invariant  sets  and  compact 
attractors . 

.  N 

If  2V  c  C([-r,0],H  )  is  a  compact  attractor  for  the 
retarded  functional  differential  equation  (RFDE) , 

(4.2)  dv(t)/dt  =  f(vfc(*)) 

then  there  is  a  Liapunov  function  V  defined  in  a  neigh¬ 
borhood  of  srf  and  satisfying  conditions  (i)-(iii)  of 
Section  2  (see  Yoshizawa  [11]). 

If  \p  €  C  (  [-r,  0]  ,  3RN  )  and  v(t,4>)  is  the  solution  of 
(4.2)  through  ip  at  t  =  0  and  $(0)(x)  =  ip(x),  -r  <_  0  £  0 
x  €  Q,  then  ip  €  C([-r,0],Xa)  and  u(t,x,^)  =  v(t,ij>), 
x  €  Q,  is  the  solution  of  (4.1)  through  ip .  Thus,  every 
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solution  of  the  RFDE  (4.2)  is  a  solution  of  (4.1). 

To  obtain  the  analogue  of  Theorem  1.1  for  (4.1),  we 

make  the  decomposition  of  Xa  as  Xa  =  U  ©  as  was 

N 

done  in  Section  2,  where  U  is  identified  with  ]R  and 

u(x)dx  . 


u  =  v  +  w,  v  €  U,  w  €  implies  v  =  |  ft  |  1 


.a 


'ft 


This  decomposition  of  X  induces  in  a  natural  way  a  de- 


composition  of  C([-r,0],X  )  as 


C([-r,0],Xa)  =  C([- r,0],]RN)  ©  C([-r,0],u£) 


<f>  =  K  +  <P,  5(0)  =  |  ft 


-1 


<j>  (9  ,x)  dx 


ft 


and 


If  u  is  a  solution  of  (4.1)  with  initial  value  u. 


u  ( t ,  x)  =  V  ( t )  +  w  ( t ,  x) 


v(t)  =  |ft|_1 

then  it  follows  from  (4.1)  that 


u (t ,x) dx 


ft 


(4.3) 

(4.4) 


V { t )  =  P(vt,wfc) 


w ( t )  =  DAw(t)  +  Q (vfc ,wfc) 


where  the  initial  values  vn,wn  v,w  are  9^ven  by 

un (0,x)dx 


0 '  0 
-1 


(4.5) 


vQ ( 9)  =  |ft|  j 


ft 


WQ (6,x)  =  uQ  (9,x)  -  v  (9) ,  -r  <  6  £  0, 


and 
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P(£,i/>)  =  | nl-1  f  f(S  +  ip(',x))dx 

Jn 

(4.6)  Q(£,*)(x)  =  f(C  +  ip(‘,x))  -  p(€,4>) 

c  e  C(  [-r,0]  ,1RN  )  ,  Il>  e  C([-r,0],U^>  . 

Equation  (4.4)  with  initial  data  can  be  written 

as 

w(t,x)  =  eDAtWg(0,x)  +  j  Q  (vg  ,wg )  ds ,  t  0 

w ( t / x)  =  wQ(t,x),  t  _<  0  , 
or,  for  -r  _<  0  £  0, 

w(t+0,x)  =  eDA^t+Q^w  (0,x)  +  [  eDA(t+0  s^q(v  w  )ds,  t+9^0, 

U  J  Q  S  S 

w(t+9,x)  =  WQ(t+0,x),  t  +  0  <  0 

If  we  define  the  n  *  n  matrix  function  X(0)  and  semi¬ 
group  S(t)  by 

Xq ( 0 )  =  0  for  0  <  0,  Xp (0)  =  I ,  the  identity 
[S ( t) Wg] (9,x)  =  eDA ^t+0^wo (0,x) ,  t  +  9  >  0 
=  wQ (t  +  0,x) ,  t  +  0  <  0 

then  the  above  formula  for  the  solution  w  becomes 

(4.7)  wt(0,x)  =  [S (t)wQ] (0,x) 

+  ft[eDA (t"s)Xn] (0)Q(vc,wc)ds 
j  0  u  s  s 
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for  all  t  ^  0,  -r£0<O.  We  will  write  this  last 
equation  as 


(4.8)  wfc  =  S (t) wQ  +  |  eD(t"s)X0Q(vs,ws)ds, 


t  >  0 


always  remembering  that  it  is  evaluated  as  in  (4.7). 

Using  the  estimates  (2.7),  equations  (4.3),  (4.8)  and 
the  same  type  of  arguments  as  in  the  proof  of  Theorem  1.1, 
one  obtains  the  following  result. 


N 

Theorem  4.1.  Suppose  Ss  <=■  C([-r,0],]R  )  is  a  compact 
attractor  for  the  RFDE  (4.2).  Then  there  is  a  6  >  0 
such  that  srf  considered  as  a  subset  of  C ( [-r ,0] ,xa) , 
a  >  3/4,  is  a  compact  attractor  for  (4.1)  if  dX  >  6 .  More 

ot 

precisely,  there  is  a  neighborhood  V  of  j/  in  C([-r,0],X  ) 
and  constants  K  >  0,  c  >  0  such  that,  for  any  Ug  6  V, 
the  solution  u(t,x,Ug),  t  >.  -r,  of  (4.1)  with 
u(0,x,uQ)  =  u0(9,x),  -r  _<  0  <  0  satisfies 


|  u ( t , • , uQ )  -  u ( t )  1  a  <  Ke-Ct ,  t  >  0 

X 


where  u (t)  =  | ft 
fies  the  equation 


i-l 


u(t,x)dx,  t  >  -r  and  u(t)  satis- 


du(t)/dt  =  f(ufc)  +  g(t,Ug),  t  >  0 


where 


I  g (t,un)  |  <  Ke-Ct , 


t  >  0. 


Examples  in  ordinary  differential  equations  are  very 
easy  to  obtain.  Conway/  Hoff  and  Smoller  [9]  have  several 
interesting  ones  for  the  case  in  which  there  is  an  in¬ 
variant  region  I .  Any  two  species  Volterra-Lotka  model 

for  which  the  ode's  have  a  unique  stable  limit  cycle  in  the 

2 

positive  quadrant  ]R+  would  have  a  compact  attractor  srf 
.  2 

in  3R+  .  Therefore,  this  model  with  large  diffusion  would 
have  the  same  attractor  srf  in  the  function  space. 

As  remarked  earlier,  it  is  almost  impossible  to  have 
an  invariant  region  when  the  equation  without  diffusion 
is  an  RFDE .  Therefore,  we  give  an  example  illustrating 
an  implication  of  the  theory  for  this  case. 

Consider  the  scalar  differential  difference  equation 

(5.1)  v ( t )  =  +  y)v(t-l) [1  +  v ( t ) ] 

in  a  neighborhood  of  y  =  0,  v  =  0.  Zero  is  always  a 
solution  of  this  equation.  Also,  for  y  <  0,  the  origin 
is  asymptotically  stable.  At  y  =  0,  the  linearized 
equation  has  two  eigenvalues  on  the  imaginary  axis  with  the 
remaining  ones  having  negative  real  parts .  For  y  >  0  there 
is  a  Hopf  bifurcation  to  a  stable  periodic  orbit.  (See, 
for  example,  Hale  [13] ,  Chow  and  Mallet-Paret  [14] ,  Stech 
[15].)  Therefore,  there  is  a  y^  >  0  and  a  neighborhood 
W  of  zero  in  C([-1,0],3R)  such  that  (5.1)  has  a  compact 


attractor  _Q^  in  W  and  ja/^  =  (0)  for  -y  ^  £  y  £0 

and  =  Wu(0)  U  y  for  0  <  y  <  yA,  where  y,,  is  the 

y  y  m  —  u  y 

periodic  orbit  obtained  from  the  Hopf  bifurcation  and 

Wu(0)  is  the  unstable  manifold  of  the  zero  solution. 

U 

Theorem  4 . 1  implies  that  is  an  attractor  for  the 

equation 

3u(t,x)/3t  =  dAu (t-x) - (^ + y ) u (t-1 ,x) [1+u (t ,x) ]  in  ft 

3u(t,x)/3n  =0  on  3ft 

in  C([-l,0,Xa),  a  >  3/4,  if  dX  >  6  where  <5  is  sufficiently 
large.  '  Furthermore,  the  estimates  in  Theorem  4.1  show,  in 
particular,  that  the  orbit  is  asymptotically  orbitally 

stable  in  C([-l,03,Xa)  if  dX  >  6. 

It  is  clear  that  the  above  remarks  remain  valid  for 
any  RFDE  in  a  neighborhood  of  a  Hopf  bifurcation.  In 
particular,  a  stable  Hopf  bifurcation  remains  stable  if 
dX  is  sufficiently  large  which  is  a  result  previously  ob¬ 
tained  by  Yoshida  [16] . 


References 


R.  Aris,  The  Mathematical  Theory  of  Diffusion  and 
Reaction  in  Permeable  Catalysts,  Clarendon  Press, 
Oxford,  1975. 

J.  Smoller,  Shock  Waves  and  Reaction-Diffusion 
Equations ,  Springer-Verlag,  1982. 

J.  Auchmuty  and  G.  Nicolis,  "Dissipative  Structures, 
Catastrophes  and  Pattern  Formation,"  Proc.  Nat. 

Acad.  Sci.,  U.S.A.  71(1974),  2748-2751. 

P.  Fife,  "Pattern  Formation  in  Reacting  and  Diffusing 
Systems,"  J.  Chem.  Phys . ,  64(1976),  554-564. 

P.  Fitzhugh,  "Mathematic  Models  of  Excitation  and 
Propagation  in  Nerves,"  Biological  Engineering, 

H.P.  Schwann,  Ed.,  McGraw-Hill,  1969. 

L.  Howard  and  N.  Kopell,  "Plane  Wave  Solutions  to 
Reaction-Diffusion  Equations,"  Studies  in  Appl . 

Math.,  52(1973),  291-328. 

R.  May,  Stability  and  Complexity  in  Model  Ecosystems, 
Princeton  Univ.  Press,  1974. 

H.G.  Othmer,  "Current  Problems  in  Pattern  Formation," 
Lecture  Notes  in  Math,  and  Life  Sci.,  Vol.  9,  1977. 

E.  Conway,  D.  Hoff  and  J.  Smoller,  "Large  Time 
Behavior  of  Solutions  of  Systems  of  Nonlinear 
Reaction-Diffusion  Equations,"  SIAM  J.  Appl.  Math., 
35(1978) ,1-16. 

D.  Henry,  Geometric  Theory  of  Semilinear  Parabolic 
Equations ,  Lecture  Notes  in  Math.,  Vol.  840,  Springer- 
Verlag,  1981. 

T.  Yoshizawa,  Stability  Theory  by  Liapunov's  Second 
Method ,  Math.  Soc.  of  Japan,  Tokyo,  1966. 

C.C.  Travis  and  G.F.  Webb,  "Existence  and  Stability 
for  Partial  Functional  Differential  Equations," 

Trans.  Am.  Math.  Soc.,  200(1974),  394-418. 

J.K.  Hale,  Theory  of  Functional  Differential  Equations 
Springer-Verlag,  1977. 


S.-N.  Chow  and  J.  Mallet-Paret ,  "Integral  Averaging 
and  Hopf  Bifurcation,"  J.  Differential  Eqns.,  26(1977) 
112-159. 

H.W.  Stech,  Hopf  Bifurcation  Calculations  for 
Functional  Differential  Equations,"  Preprint,  1983. 

K.  Yoshida,  "The  Hopf  Bifurcation  and  its  Stability 
for  Semilinear  Diffusion  Equations  with  Time  Delay 
Arising  in  Ecology,"  Hiroshima  Math.  J.,  12(1982), 
321-348. 


